Enhanced spin Hall effect by tuning antidot potential: Proposal for a spin filter 
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We propose an efficient spin filter including an antidot fabricated on semiconductor heterostructures with 
strong spin-orbit interaction. The antidot creates a tunable potential on two-dimensional electron gas in the 
heterostructures, which may be attractive as well as repulsive. Our idea is based on the enhancement of extrinsic 
spin Hall effect by resonant scattering when the attractive potential is properly tuned. Numerical studies for 
three- and four-terminal devices indicate that the efficiency of the spin filter can be more than 50% by tuning 
the potential to the resonant condition. 

PACS numbers: 72.25.Dc,71.70.Ej,73.23.-b,85.75.-d 



I. INTRODUCTION 

The injection and manipulation of electron spins in semi- 
conductors are important issues for spin-based electronics, 
"spintronics."' The spin-orbit (SO) interaction plays an im- 
portant role in the manipulation of the spins. For conduction 
electrons in direct-gap semiconductors, the SO interaction is 
expressed in the same form as that in vacuum, that is. 



so 



-cr-[px Vy(r)] 

n 



(1) 



where V(r) is an external potential and cr indicates the elec- 
tron spin s = cr/2. The coupling constant A is significantly 
enhanced by the band effect, particularly in narrow-gap semi- 
conductors such as InAs,^ compared with that in vacuum, 
\A\ - h^/{4m^f^) with mo as the electron mass and c as the 
velocity of light. 

In two-dimensional electron gas (2DEG) in semiconductor 
heterostructures, an electric field perpendicular to the 2DEG 
results in the Rashba SO interaction.^ For the electric field S 
in the z direction, the substitution of V{r) = e&z into Eq. ([T} 
yields 



^^SO = -riPyO-x- PxO-y), 

n 



(2) 



with a = eSA. Large values of a have been reported in 
experimentsji^ In the spin transistor proposed by Datta and 
. !^ ' Das^i electron spins are injected into the semiconductor het- 



erostructures from a ferromagnet, and manipulated by tuning 
the strength of Rashba SO interaction by adjusting the elec- 
tric field &. The spins are detected by another ferromagnet. It 
is well-known, however, that the efficiency of spin injection 
from a ferromagnetic metal to semiconductors is very poor, 
less than 0.1%, due to the conductivity mismatch.** To over- 
come this difficulty, the SO interaction may be useful for the 
efficient spin injection, besides the spin manipulation, in the 
spin transistor. Several spin filters have been proposed uti- 
lizing the SO interaction, e.g., three- or four-terminal devices 
based on the spin Hall effect (SHE)^'i2dLll a triple-barrier 
tunnel diode,'-' a quantum point contact,'^''^ a three-terminal 
device for the Stern-Gerlach experiment using a nonuniform 
SO interaction,'* and an open quantum dot.'^ Yamamoto and 
Kramer proposed a three-terminal spin filter with an antidot. 



using a SHE caused by the scattering of electrons at a repul- 
sive potential created by the antidot.'^ 

The SHE is one of the phenomena to create a spin current 
due to the SO interaction. There are two types of SHE. One 
is an intrinsic SHE which creates a dissipationless spin cur- 
rent in the perfect crystal. '^-^"^ Murakami et ai, for example, 
proposed that the drift motion of holes in SO-split valence 
bands induces an intrinsic SHE.'"^ The SHE of the hole sys- 
tem has been observed experimentally by Wunderlich et ai, 
using a p-n junction light-emitting diode. ^' The other type is 
an extrinsic SHE caused by the spin-dependent scattering of 
electrons by impuritiesi^2i22i2i2^ For a centrally symmetric po- 
tential around an impurity, V(r) (r - -y/x^Ty^Tz^), Eq. ([T]| 
is rewritten as 



^so 



2dV , 

-A Is, 

r dr 



(3) 



where I = (r X p)/h is the angular momentum. This results in 
the skew scattering: accompanied by the scattering from di- 
rection n to the spin is polarized in {n x n')/\n x 
In an optical experiment on the KeiT rotation, Kato et al. ob- 
served a spin accumulation at sample edges along the electric 
current in «-type GaAs^ which is ascribable to the extrin- 
sic SHE. The experimental result has been quantitatively ex- 
plained by a semi-classical theory considering the skew scat- 
tering and "side jump" effects. 

In our previous paper,^^ we have quantum-mechanic ally 
formulated the extrinsic SHE for 2DEG in semiconductor het- 
erostructures. We have examined the SHE due to the scatter- 
ing by an artificial potential created by an antidot, scanning 
tunnel microscope (STM) tip, etc. An antidot is a small metal- 
lic electrode fabricated above 2DEG, which creates an elec- 
trically tunable potential on 2DEG. The potential is usually 
repulsive, but it could be attractive when a positive voltage is 
applied to the antidot. We have found that the SHE is sig- 
nificantly enhanced by resonant scattering when the attractive 
potential is properly tuned. 

We have stressed that the extrinsic SHE is easier to under- 
stand in 2DEG than in three-dimensional c ase. Let us consider 
an axially symmetric potential V(r) (r - -¥ y^) created by 
an antidot on conduction electrons in the xy plane. The SO 
interaction in Eq. ([U becomes 



Hso^-A-'^ls, = Vi{r)ks„ 
r dr 



(4) 
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FIG. 1: Our model for (A) three- and (B) four-terminal devices for 
the spin filter. They are fabricated on two-dimensional electron gas 
in the xy plane. Both the devices include an antidot at the center of 
junction, which is a square area surrounded by broken lines. Three 
or four ideal leads connect the junction to reservoirs. Reservoir 1 
is a source from which unpolarized electrons are injected into the 
junction. The voltages are equal in the drains; reservoirs 2 and 3 in 
model (A) and reservoirs 2, 3, and 4 in model (B). 



where L and are the z component of angular momentum and 
spin operators, respectively. Vi(r) - -(2A/r)V'(r), which has 
the same sign as V{r) if \V(r)\ is a monotonically decreasing 
function of r and A > 0. Assuming that V(r) is smooth on 
the scale of the lattice constant, we adopt the effective mass 
equation 



-- — A + V(r) + Vi(r)y. 
2m* 



(5) 



for an envelope function if/{r) with m* as the effective mass. 
In Eq. (|5]l, and s~ are conserved in contrast to the three- 
dimensional case with Eq. (O. For 5- - +1/2, an electron 
experiences the potential of V{r) + Vi(r)l~/2. As a conse- 
quence, the scattering of components of I- > Q (L < 0) is 
enhanced (suppressed) by the SO interaction for = 1/2 
when Vi{r) has the same sign as V{r). The effect is oppo- 
site to that for s. - -1/2. This is the origin of the extrinsic 
SHE. We have formulated the SHE in terms of phase shifts in 
the partial wave expansion for 2DEG and shown that the SHE 
is largely enhanced by resonant scatteringi^S These results are 
summarized in Appendix A. 

In the present paper, we consider three- and four-terminal 
devices including an antidot, as shown in Fig.lT] and examine 
the enhancement of the SHE. We evaluate an efficiency of the 
spin-filtering effect by resonant scattering in the case of at- 
tractive potential. Although our three-terminal device is very 
similar to the spin filter proposed by Yamamoto and Kramer, 
they have only studied the case of repulsive potential. We 
show that our device can be a spin filter with an efficiency of 
more than 50% by tuning the potential to the resonant condi- 
tion. 

We numerically solve the effective mass equation [Eq. ©I 
with an appropriate boundary condition for our devices. A 
confining potential for the leads (quantum wires) could induce 
the SO interaction, following Eq. ([X]| 30,3i,32.3 3 jjowever, its ef- 
fect on the electrons is much smaller than the SO interaction 
induced by the antidot potential because the amplitude of the 
wavefunction is small around the edges of the leads. There- 



fore, we consider the antidot-induced SO interaction only. We 
also assume that the antidot potential V(r) is independent of 
z. Otherwise, it would create the Rashba-type SO interaction, 
Eq. dill with a = A(dV/dz), in addition to Eq. ©. The Dres- 
seUiaus SO interaction is also disregarded, which is induced 
by the inversion asymmetry of the crystali^ These effects will 
be discussed in the final section. 

The electron-electron interaction is not taken into account. 
The Coulomb blockade is not relevant to the case of antidot, 
in contrast to that of conventional quantum dot which is con- 
nected to the leads via tunnel barriers.^'' In our model, there- 
fore, the influence of the electron-electron interaction is only 
quantitative and could be considered by the mean-field level, 
as discussed in the final section. Note that there have been 
several researches for the spin-cuiTent generation based on the 
electron-electron interaction in the absence of SO interaction, 

e.g., using single or double quantum dot a^^'^^i^^ and quantum 
wires.3''-4Q.4i,42,43 

The organization of the present paper is as follows. In Sec. 

II, we describe our model for three- and four-terminal devices 
and calculation method. The calculation of spin-dependent 
conductance in multi-terminal devices is formulated using the 
Green's function in the tight-binding model. ii^^i^^i^ In Sec. 

III, we present numerical results of the conductance and spin- 
filtering effect when the strength of antidot potential is tuned. 
We also investigate the density of states (DOS) in the junction 
area of the devices to illustrate vutual bound states. The exis- 
tence of the virtual bound states at the Fermi level is a strong 
evidence that a resonant scattering takes place when the spin- 
filtering effect is remarkably enhanced. In addition, we per- 
form a channel analysis of the spin-dependent conductance to 
closely examine the resonance. The final section (Sec. IV) is 
devoted to the conclusions and discussion. 



II. MODEL AND CALCULATION METHOD 

In this section, we explain our model and calculation 
method. We numerically solve the effective mass equation in 
the tight-binding model for the devices. In the presence of the 
SO interaction in Eq. (|4|i, the z component of electron spin s. 
is conserved although l~ is not a good quantum number owing 
to the lack of rotational symmetry in our devices. Hence we 
can solve the equation for = +1/2 separately. 



A. Model 

We consider three- and four-terminal devices with an anti- 
dot, fabricated on semiconductor heterostructures, as shown 
in Fig.[T] Three or four leads (quantum wires) of width W are 
joined to one another at a junction, which is a square area sur- 
rounded by broken lines in the figure. The leads are formed by 
hard-wall potential and connected to the reservoirs. Reservoir 
1 is a source from which unpolarized electrons are injected 
into the junction. The electrons are outgoing to the drains; 
reservoirs 2 and 3 (2, 3, and 4) in the three-terminal (four- 
terminal) device. The voltages are equal in all the drains. 
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An antidot creates an axially symmetric potential V{r), 
where r is the distance from the center of the junction. It is 
assumed to be attractive and given by a smooth potential well, 



V(r) = 



Vo 



^{l-sin(;r^)} (\r - Ro\ < ^) (6) 







ir-Ro>^) 



with Vo < 0. The radius of the potential well is Ro - W/4, and 
we choose ARq - O.IRq. The gradient of V{r) gives rise to the 
SO interaction in Eq. Q. 

For the numerical study, we discretize the two-dimensional 
space with lattice constant a (tight-binding model with a 
square lattice). The width of the leads is W = {N + l)a with 

= 29; the wavefunction becomes zero at the zeroth and 
(A^ + l)th sites in the transverse direction of the leads. The 
Hamiltonian is 



H 



+ i,i;i,i+\;cr'^ i,j,a'' i.j+\,a- + c), 



(7) 



where and c are creation and annihilation operators 

i.j\(T i,j;cr f 

of an electron at site (/, j) with spin cr. Here, t = h^/2m*a^, 
where m* is the effective mass of electrons. Vjj is the potential 
energy at site j), in units of t. The transfer term in the x 
direction is given by 

Tij-i+ij-± - 1 + iA{Vi+i/2j+i - Vi+i/2j-i), (8) 
whereas that in the y direction is 

Tij;ij+i-± - 1 + iMVi+ij+i/2 - Vi-ij+i/i), (9) 

where A - A/ (4a^) is the dimensionless coupling constant of 
the SO interaction. V,+i/2,j is the average of the potential en- 
ergy at sites (;, j) and (; + l,j), and V,-,)+i/2 is that at sites (/, j) 
and (ij+l). 

The SO interaction is absent in the leads. The wavefunction 
of conduction channel yu in the leads is written as 



Ufi(i') 



exp{ik^aj')u^(i'), 
njii' 



N+ 1 



sm 



N+ 1 



(10) 
(11) 



where /' and / denote the site numbers in the transverse 
and longitudinal directions of the leads, respectively. The 
wavenumber k/j is determined from the condition that the dis- 
persion relation 

Ef,ik) = 4t- 2f cos ^ ^"^ ^ j - 2tcos(ka) (12) 

is identical to the Fermi energy Ep. The band edge of channel 
fi is defined by Efj(k - 0). The band edge is located below Ep 
for the conduction channel. For E^{k = 0) > Ep, on the other 



hand, mode fi is an evanescent mode. The wavefunction of the 
evanescent mode is given by 



^iiii'J') = exp(-Af^fl/)M^(/'), 



(13) 



where a f is the distance from the junction along the lead and 
satisfies E^iiK^^) - Ep. 



B. Calculation method 

For numerical studies, we introduce the Green's function. 
The Green's function for the junction area is defined by 

-1 



CAE) = 



EI-'H^-J^ I." 



(14) 



where Wo- is the truncated matrix of the Hamiltonian for the 
junction area (A^^ x A^^) with spin cr, and I.'' is the self-energy 
due to the coupling to the lead p: 



Here f/ is a unitary matrix, U - (mi, M2 
[m^(1),m^(2),--- ,M^(A^)]'. A = diag(.li,.l2 



(15) 

Un), where = 
, ■ ■ ■ , Aff), where 
Ajj - exp(ikf,a) for conduction channels and Ajj - exp(-/f^a) 
for evanescent modes. Tp is a coupling matrix (A^ x A^^) be- 
tween lead p and the junction: Tp{pi, i) = 1 for being an 
adjacent site in lead p to site / in the junction; Tp{pi, i) = 
otherwise.'*'* 

The spin-dependent conductance from reservoir p to reser- 
voir q is obtained from the Landauer-Biittiker formula. It is 
written as 

2 

G^P = ^Tr [r'GAEW'GUE)] , (16) 



where 



h 



TP ^i[I.P 



The total conductance is CP = G'f + G'"' . 
polarization in the z direction is defined by 

f^qp ^qp 

CP ^ W -G_ 



(17) 

whereas the spin 
(18) 



for the current from reservoir p to q. 

To elucidate the virtual bound states in the potential well, 
we calculate the DOS in the junction area. It is evaluated from 
the Green's function (fT4l i as"*^ 



D{E) ^ --Y^lm\YxGAE)\. 



(19) 



We assume that A - Q.l for the strength of SO interaction, 
which corresponds to the value for InAs, A - 117.1 A-,- with 
a — WjliQ) and width of the leads W ~ 5Q nm. The temperature 
r = 0. We focus on the transport from reservoir 1 to 2 and 
omit the superscript 21 of G^' and P?', otherwise stated. Note 
that the conductance from reservoir 1 to 3 is related to G|' = 
G^' from the symmetry of the system, for both three- and four- 
terminal devices. The current from reservoir 1 to 4 is not spin- 
polarized in the four-terminal device. 
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FIG. 2: (Color online) Numerical results for the three-terminal de- 
vice with kpRo = 2, where Ro is the radius of attractive potential. (A) 
Conductance G± from reservoir 1 to 2 in Fig.mA) for s, = ±1/2 
and (B) spin polarization P- of the output current in reservoir 2, as 
functions of the potential depth |Vo|. In (A), solid and broken lines 
indicate G+ and G , respectively. A dotted line shows the conduc- 
tance per spin in the absence of the SO interaction. (C) Grayscale 
plot of the density of states in the junction area, D(E), in the plane of 
I Vol and energy E of electron. 



III. CALCULATED RESULTS 



We calculate the conductance G± for spin s. - ±1/2 and 
spin polarization P, when the potential depth IVqI is tuned. 
We examine three cases of ^p^o = 1, 2, and 3, where the 
Fermi wavenumber kp is defined by the Fermi energy Ep as 
Ep/t - (kpa)^.*'^ In the three cases, the Fermi energy Ep is 
different, whereas the radius of the potential well Rq is fixed. 
The number of conduction channels in the leads is 1,2, and 3, 
respectively. 

Here we discuss the cases of kpRo = 2 and 3. The numerical 
result with kpRo = 1 is given in Appendix B. (Surprisingly, we 
find a perfect spin polarization P. = 1 in the case of kpRi) - 1. 
However, the transport properties seem quite specific. This is 
due to a strong interference effect in the case of single con- 
duction channel.) 



A. Caseof ytp/fo = 2 

We begin with the three-terminal device in the presence of 
two conduction channels in the leads (kpRo = 2). Figures|2jA) 
and|2]!B) show the conductance G± for 5, = ±1/2 and spin 
polarization P., respectively, when the potential depth IVqI is 
gradually changed. As seen in Fig.|2jA), the conductance G± 
shows three minima as a function of IVqI. At the first minimum 
at \Vo\/Ep Si 0.6, the difference in the conductance for = 
±1/2 is small. At the second and third minima at \ V{)\/Ep ^ 2 
and 5, respectively, the difference is remarkable, which results 
in a large spin polarization in the z direction [Fig. |2tB)]. P^ 
is enhanced to 25% around the second minimum and 61% 
around the third minimum. 

The behavior of G± should be ascribable to resonant scat- 
tering at the potential well. The resonant scattering takes place 
through a virtual bound state in the potential well, which en- 
hances the electron scattering to the unitary limit (Appendix 
A). This makes the minima of G± in our situation of both 
three- and four-terminal devices. (It is not trivial whether the 
resonant scattering makes minimum or maximum of the con- 
ductance in multi-terminal devices. See the discussion in Ap- 
pendix A.) Around the minima of the conductance, the differ- 
ence between G+ and G is magnified. Look at the third min- 
imum of the conductance around \Vo\/Ep ^ 5. The minimum 
of G+ is located at a smaller value of |Vo| than that of G . In 
consequence, P^ shows a pair of negative dip (G+ < G_) and 
a positive peak (G+ > G_). This dip-peak structure of P, can 
be understood in terms of phase shifts when the resonance for 
s~ - +1/2 is well separated from each other (Appendix A). 

To confirm the above-mentioned scenario regarding the res- 
onant scattering, we calculate the DOS in the junction area. 
Figure HJC) shows a grayscale plot of the DOS D{E) in the 
plane of IVqI and energy E of electron. The band edge for the 
lowest and second conduction channels in the leads [Ei(k = 0) 
and E2(k = 0) in Eq. ^] are located at E/Ep = 0.154 and 
0.615, respectively. The sharp peaks of D(E) below the low- 
est band edge coiTespond to the bound states inside the junc- 
tion area. With an increase in the potential depth IVol, several 
bound states appear one after another. The first one is an 5 - 
like bound state (I. = 0) although the angular momentum I. is 
not a good quantum number in our device because of the lack 
of rotational symmetry. The bound state exists even without 
the potential well (IVol - 0) in the junction area,"^ and changes 
to the S -like bound state in the potential well with increasing 
IVol. The state is doubly degenerate due to the Kramers degen- 
eracy. The next are P-like bound states (L = ±1). They are a 
pair of Kramers degenerate states. Then D-like bound states 
(L = ±2) appear, which are clearly split into two by the SO 
interaction. Another S -like state is located at approximately 
the same energy. Finally F-liks bound states (L - ±3) appear 
in Fig.|2|C). The pair of Kramers degenerate states is largely 
separated for the F-liks states. 

The peaks of the bound states in D(E) are broadened above 
the band edge of the lowest conduction channel in the leads, 
which significantly influence the electron scattering at the 
Fermi level as virtual bound states. The second and third min- 
ima of the conductance G+ are located at the position of D 
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FIG. 3: (Color online) Numerical results for the four-terminal device 
with kj:R(, = 2, where Rq is the radius of attractive potential. (A) 
Conductance G± from reservoir 1 to 2 in Fig.[nB) for j; = ± 1 /2, (B) 
conductance G+' from reservoir 1 to 4, and (C) spin polarization 
of the output current in reservoir 2, as functions of the potential depth 
|Vo|. In (A) and (B), solid and broken lines indicate G+ and G_, or 
G\' and G1' , respectively. A dotted line shows the conductance per 
spin in the absence of the SO interaction. In (B), solid and broken 
lines for G+' are completely overlapped. (D) Grayscale plot of the 
density of states in the junction area, D(E), in the plane of |Vo| and 
energy E of electron. 



and F-like virtual bound states at Ep, respectively. This is a 
clear evidence of the resonant scattering through virtual bound 
states. (At the first minimum of G± around \Vo\/Ep = 0.6, we 
cannot find any virtual bound state at the Fermi level. The 
minimum of G± may not be due to the resonant scattering, but 
due to some interference effect around the junction.) 

We present the calculated results for the four-terminal de- 
vice with ^F^o = 2 in Fig. [3j (A) conductance G± for 
s, = +1/2 from reservoir 1 to 2 in Fig. [HB); (B) conduc- 
tance G^' from reservoir 1 to 4, and (C) spin polarization P- 
of the output current in reservoir 2, as functions of the poten- 
tial depth IVol. As seen in Fig.[3lB), G^' = Gl^ because the 
SHE does not make a spin polarization in the output current in 
reservoir 4. The characters of conductance G± for s. - +1/2 
and spin polarization P, are almost the same as those in Fig.|2] 
for three-terminal device. The conductance shows three min- 




FIG. 4: (Color online) Numerical results for the three-terminal de- 
vice with kpRg = 3, where Rg is the radius of attractive potential. (A) 
Conductance G+ from reservoir 1 to 2 in Fig. [HA) for = ±1/2 
and (B) spin polarization P- of the output current in reservoir 2, as 
functions of the potential depth | Vol- In (A), solid and broken lines in- 
dicate G+ and G_, respectively. A dotted line shows the conductance 
per spin in the absence of the SO interaction. (C) Grayscale plot of 
the density of states in the junction area, D(E), in the plane of | Vol and 
energy E of electron. Regarding the result for four-terminal device 
with kpRg = 3, a broken line in (B) indicates the spin polarization P- 
of the output current in reservoir 2 in Fig.[TjB). 



ima. The second and third minima are clearly due to resonant 
scattering via D- or F-like virtual bound states, as seen in the 
DOS in Fig. [3jD). Around the minima, the conductance for 
= +1/2 is largely split by the SO interaction, which results 
in a large spin polarization P^. 



B. Caseof ^f/^o = 3 

Figure |4] shows the calculated results for the three-terminal 
device in the presence of three conduction channels in the 
leads {kpRo = 3): (A) conductance G± for = ±1/2 and 
(B) spin polarization P^ in the z direction, as functions of the 
potential depth IVqI. Figure|4|C) shows a grayscale plot of the 
density of states D(E) in the plane of IVqI and energy E of elec- 
tron. In Fig.HtB), a broken line indicates the spin polarization 
in the four-terminal device. 

The conductance G± shows several minima as a function 
of potential depth IVqI- The spin polarization P^ is enhanced 
around the minima of G±. These properties can be understood 
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FIG. 5: A channel analysis for incident waves from reservoir 1 in the 
three-terminal device with kpRg = 2. (A) Spin polarization P- of the 
output current in reservoir 2 in Fig.[TjA), as a function of the poten- 
tial depth I Vol, for the incident electrons in the lowest channel (solid 
line) and second channel (broken line). (B) Conductance G+ from 
reservoir 1 to 2 for the incident electrons in the lowest channel with 
= 1/2 (solid line) and -1/2 (broken line). A dotted line indicates 
the conductance per spin in the absence of the SO interaction. 



in the same way as in the preceding subsection. The polar- 
ization p. is enhanced to 41% at \Vo\/Ep - 3.1 in the three- 
terminal device, and it is enhanced to 49% at \Vo\/Ep = 3.2 in 
the four-terminal device. This is due to resonant scattering via 
G-like virtual bound states {L - ±4). 

Compared with the case of two conduction channels in 
Figs. |2] and |3] the values of the conductance G± are larger in 
the case of three conduction channels {kpRo - 3), whereas the 
maximum value of spin polarization is almost the same. This 
implies a more efficient spin filter in the case of three conduc- 
tion channels than in the case of two conduction channels. 



C. Channel analysis for spin filtering 

In cases of kpRo - 2 and 3, there are two and three con- 
duction channels in the leads, respectively. To examine the 
resonant scattering in detail, we perform a channel analysis of 
incident waves from reservoir 1 . The results are given only 
for the three-terminal device in this subsection. 

In the case of kpRo = 2, we plot the spin polarization 
for the incident electrons in the lowest and second channels in 
Fig.|5jA). At IVol/Ef ~ 2 (resonance by D-like virtual bound 
state), P, is enhanced to 73% for the lowest channel while 
it is to 18% for the second channel. Hence the former plays 
a main role in the spin polarization. At \Vo\/E-p ~ 5 (reso- 
nance by F-like virtual bound state), on the other hand, 
becomes 75% for the lowest channel while it becomes 83% 
for the second channel. Then both channels are important for 
the spin-dependent scattering. 

We could selectively inject the lowest channel to the junc- 
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FIG. 6: A channel analysis for incident waves from reservoir 1 in 
the three-terminal device with kpRo = 3. (A) Spin polarization P; 
of the output current in reservoir 2 in Fig.[TjA), as a function of the 
potential depth \ Vo\, for the incident electrons in the lowest channel 
(solid line), second channel (broken line), and third channel (dotted 
line). (B) Conductance G+ from reservoir 1 to 2 for the incident 
electrons in the lowest channel with s- = 1/2 (solid line) and -1/2 
(broken line). A dotted line indicates the conductance per spin in the 
absence of the SO interaction. 



tion, e.g., using a quantum point contact fabricated on the lead 
connected to reservoir 1 . Then we could realize a spin filter 
with an efficiency of about 75%. In Fig.|5jB), we plot the con- 
ductance G± when only the lowest channel is injected from 
reservoir 1. At \V[)\/Ep ~ 2, the conductance almost vanishes 
although P; is enhanced to 73%. This is due to an interfer- 
ence effect at the junction as in the case of single conduction 
channel with kpRo = 1 (Appendix B). At \ Vo\/Ep ~ 5, on the 
other hand, the total conductance is G+ + G- = Q.4{e^/h) and 
= -75%. The latter situation is favorable to application to 
a spin filter 

A similar channel analysis is given for the case of ^f^o = 3 
in Fig. |6l There are three incident channels in this case. It 
is notable that, at \V[)\/Ep ~ 2.8, a spin polarization of P, = 
62% is realized for the incident electrons in the lowest channel 
while the total conductance is as large G_ = OMe^/h). 



D. Repulsive potential 

We investigate the SHE caused by the scattering by a repul- 
sive potential, Vq > in Eq. Figure |7] shows (A) conduc- 
tance G± for i, - +1/2 and (B) spin polarization P^ in the z 
direction, when the potential height Vq is gradually increased. 
The extrinsic SHE is expected even with a repulsive poten- 
tial.'^ However, the spin-filtering effect is much weaker than 
the case with an attractive potential. In Fig. IllB), the spin 
polarization is at most P, x 0.3% in the three-terminal de- 
vice and a: 0.45% in the four-terminal device. In this case, 
the resonant scattering does not take place since virtual bound 
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FIG. 7: Numerical results for the three-terminal device (curves la- 
beled by a) and four-terminal device (curves labeled by b) with re- 
pulsive potential [Ko > in Eq. kpRo = 2, where Rq is the radius 
of the potential. (A) Conductance G± from reservoir 1 to 2 in Fig. 
[T] as a function of potential height Vo, for = 1/2 (solid lines) and 
-1/2 (broken lines). (B) Spin polarization P- of the output current in 
reservoir 2. 

states are hardly formed in the potential barrier. This indicates 
an important role of resonant scattering in the enhancement of 
the SHE discussed in the preceding subsections. 



IV. CONCLUSIONS AND DISCUSSION 

We have numerically studied the extrinsic SHE in multi- 
terminal devices including an antidot, fabricated on semicon- 
ductor heterostructures with strong SO interaction. The an- 
tidot creates a tunable potential on two-dimensional electron 
gas in the heterostructures, which may be attractive as well 
as repulsive. When an attractive potential is tuned properly, 
the resonant scattering via a virtual bound state takes place, 
which makes minima of the conductance from reservoir 1 to 
2 in Fig.[I] Then the difference between the conductances for 
s, = ±1/2 is enlarged and, as a result, the spin polarization 
is significantly enhanced in the direction perpendicular to the 
two-dimensional plane. The spin polarization can be more 
than 50% in our three- and four-terminal devices. 

The enhancement of the extrinsic SHE by resonant scatter- 
ing has been studied in different systems. Kiselev and Kim 
have proposed a three-terminal spin filter without antidot in 
the presence of Rashba SO interaction [Eq. (|2]i].'°'" They 
have pointed out that the spin-filtering effect is enhanced by 
resonant scattering at the junction area when the Fermi energy 
of 2DEG is tuned. In their device, the direction of spin polar- 
ization is tilted from the z direction perpendicular to the plane. 
In our device, the spin is polarized in the z direction, which is 
easier to detect by an optical experiment on the Kerr rotation^** 
and, above all, more suitable to spintronic devices. 

The extrinsic SHE enhanced by (many -body) resonant scat- 



tering has been examined for metallic systems with magnetic 
impurities.i^iJiS In the case of semiconductor heterostruc- 
tures, however, we have a great advantage in the tunability of 
scattering potential. The enhanced SHE caused by resonant 
scattering at a single potential can be investigated in details. 

We make some comments regarding our calculations, (i) 
The electron-electron interaction has been neglected. Below 
the band edge of the lowest conduction channel in the leads, 
we have observed bound states in the density of states in the 
potential well. Since the bound states are occupied by elec- 
trons, we have to consider the electron-electron interaction 
between the electrons and conduction electrons at the Fermi 
level. The Hartree potential from the trapped electrons should 
be taken into account although the Coulomb blockade is irrel- 
evant to the case of antidot potential without tunnel barriers, 
in contrast to the case of conventional quantum dots.^^ In our 
calculated results, therefore, the values of | Vol at the resonance 
are underestimated. 

(ii) It is necessary to create such a deep potential as IVol ~ 
Ep in designing the devices. This might be difficult with 
a usual antidot structure fabricated on semiconductor het- 
erostructures. Alternatively, we could make such a deep po- 
tential using a STM tip, a charged impurity under the antidot, 
etc. 

(iii) We have assumed that the antidot potential V{r) is inde- 
pendent of z- Otherwise, the Rashba-type SO interaction, Eq. 
^ with a = A{dV/dz), must be added to Eq. (|4]l. This would 
create an effective magnetic field in the xy plane and thus de- 
crease the spin polarization in the z direction. The Dresselhaus 
SO interaction has also been disregarded. The SO interaction 
is induced by the inversion asymmetry of the crystal^ and ex- 
pressed as 

^^so = ^i-p.xcr., + PyO-y). (20) 

This would also result in an effective magnetic field in the xy 
plane and lessen the spin polarization in the z direction. 
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APPENDIX A: FORMULATION OF SPIN HALL EFFECT IN 

2DEG 

Here, we summarize our previous study in Ref . 

First, we give a quantum mechanical formulation of the ex- 
trinsic SHE for 2DEG. For the scattering problem with Eq. 

we adopt the partial wave expansion with L - m - 0, 
±1, ±2, ■ ■ ■ r'^ As an incident wave, we consider a plane wave 
propagating in the x direction, e'*^*, with spin - 1 12 or 



-Ijl. E - H k /(2m*). The plane wave is expanded as 



^ikx _ ^ikr cos 6 



2 rJmikr)e'"'', (Al) 



where 6 is the angle from the x direction and is the mth 
Bessel function. Its asymptotic form at r — > oo is given by 
Jm{kr) ~ V2/ (nkr) cos(kr - mn/l - 7t/4). In the solution of 
Eq. (|5j, Jm(kr) in Eq. ( lAll i is replaced by for s, - ±1/2, 



(A2) 



for spin = +1/2, where /J^(r) satisfies 



2m* \ c/r^ r dr i '' 



= ERlir). (A3) 
Its asymptotic form determines the phase shift 5* : 

^'"^'"^ ~ VJ;^"' - T - ^ ^ ■ ^^^^ 

From Eqs. ( IA3l l and (IA4l i. we immediately obtain the relation 
of = indicating the time reversal symmetry. The SO 
interaction is not relevant to the S wave {m - 0): 6^ = <5q s 

The scattering amplitude f^{6) for spin s, - ±1/2 is ex- 
pressed in terms of phase shifts. The asymptotic form of the 
wavefunction in Eq. (IA2l i is given by 



J[kr+nlA) 



where f'^iff) is related to the S matrix by 
tie) = Y.f^ne""', 

m 

From these equations, we obtain 

/^(0)=/l(0)±/2(0), 



(A5) 

(A6) 
(A7) 

(A8) 



where 



i YZnk 

oo 



m=l 



fiie) = 



1 °° 



(A9) 
(AlO) 



The scattering cross section is given by cr'^iO) - 
Hence the spin polarization of the scattered wave in the 6 di- 
rection is expressed as 



ri' - ^ 2Re(/i/2*) 

I/+P + I/-P I/1P + I/2P' 



(All) 
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FIG. 8: Partial wave expansion for the extrinsic SHE due to the scat- 
tering by a potential well, Eq. ( IA13t . in 2DEG. kR^ = 1 . The strength 
of the SO interaction is Ak- = 0.01. (A) Scattering probability of 
each partial wave, sin^ dj;„ (B) scattering cross section cr*(6 = -n/2) 
for .S; = ±1/2, and (C) spin polarization al 9 = -n/2, as 
functions of the potential depth |Vol [normalized by electron energy 
E = Ji^k^ /(2m')]. In (A), solid and broken lines indicate the cases 
of s, = 1/2 and -1/2, respectively, for m = 1 ((5*| = i5p. The 
scattering probability for \m\ > 2 is negligible. In (B), a broken line 
indicates the cross section aid = -n/2 in the absence of SO interac- 
tion. Inset: scattering cross section in the absence of SO interaction 
at = -0A51T (a), -n/2 (broken line), and 6 = -0.55n (b). 



when the incident electron is unpolarized. This formula is 
analogous to that of skew scattering in three dimensionsi^^ 
The spin is polarized in the z direction and Pri-ff) - -P,{ff). 

A remark is made on Eq. ( IA5b . We put a phase of n/A on 
the exponent of the second term on the right side. Although 
the phase has no physical meaning, it ensures the "optical the- 
orem" that the total cross section is determined only by the 
amplitude of forward scattering; 



pin 



-Im/^(0). 



(A12) 



Now we apply our formula of the SHE to the case of a po- 
tential well, Eq. ^ with A/^o -> 0, that is. 



V{r) = Va9(RQ - r) 



(A13) 



(Vo < 0), where Oit) is a step function [9(1) = 1 for f > 0, for 
f < 0]. Then V, = -(2A/r)V'(r) = (2A/Ro)Vo6{r - Rq) with 
the 6 function 6{t). 

The phase shifts 6% are calculated by solving Eq. iA3i . The 
calculation is elementary. We have only to consider the case 
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FIG. 9: Partial wave expansion for tlie extrinsic SHE due to the scat- 
tering by a potential well, Eq. l lA13t , in 2DEG. kRg = 2. The strength 
of the SO interaction is Ak' = 0.01. (A) sScattering probability of 
each partial wave, sin^ (5*, and (B) spin polarization P-at6 = -n/2, 
as functions of the potential depth |Vol [normalized by electron en- 
ergy E = h^k^ /(2m')]. In (A), solid and broken lines indicate the 
cases of =1/2 and -1/2, respectively, for ra > ((5*„, = Sf„). 



with m > because of the relation of S%„ = 6"^^. For r > R^), 
V - and thus 



R%(r) = CiJ„,(kr) + C2Y,„(kr) 
y nkr' 

+C2 &in(kr - mnjl - 7r/4)], 



[Ci cos(A;r - mnl2 - k/4) 



(A14) 



(A15) 



where y,„ is the mth Neumann function. From Eqs. ( IA15b and 
(lA4l i. we find 



tandt, ^ -C2/C1. 
For r<Ro,V^ Vq. Then 



(A16) 



(A17) 
(A18) 



cross section cr*(9 = -7r/2) for s- - +1/2, and (C) spin po- 
larization P, at = -?r/2, when the potential depth IVqI is 
gradually changed. The strength of the SO interaction is set 
to be Ak^ - 0.01, which corresponds to the value for InAs, 
/I = 117.1 A^,^ with electron wavelength 27r/A; - 70 nm. The 
radius of the potential well is Rq - 1/k. As seen in Fig. [St A), 
with an increase in IVqI, the scattering probability increases 
and becomes unity at some values of | Vol (unitary limit with 
6% = 7r/2) for m - (S wave) and m - +1 (P wave). This 
is due to resonant scattering through virtual bound states in 
the potential well. The resonant width is narrower for larger 
\m\ because of the centrifugal potential oc m^/r^ separating the 
bound states from the outer region. 

As discussed in the text with Eq. (|5]l, the scattering of a 
partial wave with positive m is enhanced (suppressed) for spin 

= 1/2 (ij = -1/2). Around the resonance of the P waves, 
the scattering of (m, Sj.) = ( 1 , 1 /2) goes to the unitary limit at a 
smaller value of IVqI than the scattering of (m, s^) - (1,-1/2). 
Thus we observe a difference between 5| and in Fig.[8jA). 
This leads to the difference of the scattering cross section at 
6 = -7r/2 for = +1/2 [Fig.[8jB)] and a spin polarization 
[P, X 30% in Fig. HO]. When 6^ {\m\ > 2) is negligible, 
Eqs. (lASTl- dATOl l yield 



0-^(0 = -7T/2) = — -[sin^ 6q + sin^ A^i 
nk 

±2 sin^o sin(25i - 60) sin A5i], 



(A21) 



where A61 — 6'^ - 6^ and 61 - (5| + 6\)I2. Around the 
resonance, A^i is enlarged, which results in the enhancement 
of the SHE. For 6\ ~ 7r/2, the spin polarization is given by 



P,{e = -7r/2) 



2 sin" 6q sin A5i 
sin^ 6q + sin^ A^i 



(A22) 



from Eq. (lAlll i. 

In our three- and four-terminal devices, the resonant scat- 
tering makes minima of the conductance. It should be hard 
to say, however, whether the resonant scattering enhances or 
suppresses the conductance in general. For simplicity, let us 
ignore the SO interaction. We plot the scattering cross section 
in the inset in Fig.[8jB), in the direction ofO — -0.457r, -7r/2, 
and 9 - -O.SSn. Using 60 and 6^=6^=61, the cross section 
is written as 



The connection of Eqs. ( IA14l i and ( IA17l i at r = Rq yields 
[J„-i(kRQ) - J„+i(kRQ)]JJk'Ro) - aMkRo) 



tan^„ 



[Y^^iikRo) - Y„+i(kRo)]J„,{k'RQ) - at,Y„ma) 

(A19) 



with 



ai = {k'lk)[J„^i{k'Rt;) - J„^i{k'Ro)] 
+2m[l + (k'/kf](kA/Ro)J„,(k'Ro). 



(A20) 



Here, we have used the relation of dJ,„(x)/dx - [7,„_i(x) - 
Jm+i(x)]/2 and dY„,(x)/dx = [Y„-i(x) - Y„+i(x)]/2. 

The calculated results are shown in Fig. |8] (A) scatter- 
ing probability of each partial wave, sin^^,"^,, (B) scattering 



o-(6) - — [sin^ So + 4 sin^ 61 cos^ 
nk 

+4 sin^osin^i cos((5o - 5i)cos0]. 



(A23) 



The resonance of the P wave (61 - 7t/2) has no effect on the 
scattering at = -7r/2. Its effect on the scattering at 6 ^ 
-7r/2 depends on the value of Sq. In the inset in Fig.OB), the 
resonance of the P wave makes a peak of cr{6 = -0.457r) and 
a dip of cr(6 = -0.55;7r). In the situation of our devices, the 
conductance should reflect an integrated value of cr(0) around 
9 - -7r/2, and also a complicated interference effect. 

Figure|9]shows the calculated results in the case of kR^ - 2. 
The resonant scattering is observed for < \m\ < 3. Around 
the resonance of F waves (|m| - 3), P, is enhanced to 72%. 
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where Adj, = 6'^-6^ and ^3 = (6^ 
(3, 1/2), 6 



FIG. 10: (Color online) Numerical results for the three-terminal de- 
vice with kjiRd = 1, where Rq is the radius of attractive potential. (A) 
Conductance G+ from reservoir 1 to 2 in Fig. [TjA) for s- = ±1/2 
and (B) spin polarization P, of the output current in reservoir 2, as 
functions of the potential depth |Vo|. In (A), solid and broken lines 
indicate G+ and G_, respectively. A dotted line shows the conduc- 
tance per spin in the absence of the SO interaction. (C) Grayscale 
plot of the density of states in the junction area, D(E), in the plane of 
I Vol and energy E of electron. 



In general, a sharper resonance enlarges 6^ - 6„ for larger \m\, 
which results in a larger polarization. 

In Fig. IUB), P- shows a dip and peak structure at the F- 
wave resonance. This structure is observed when the reso- 
nance of (m, sS) - (3, 1/2) is sufficiently separated from that 
of (;«, s,) = (3,-1/2). Around the F-wave resonance, (5| ~ 
and 62 ~ 62 (= ^2), as seen in Fig. IUA). Neglecting ^0, Eq. 



dAllb yields 



-6^)/2. At the resonance of 
(m, s,) = (3, 1/2), 5+ = 7r/2 and ^3 ^ 0. Then P,(9 = -n/2) 
shows a dip for 7r/2 < 62 < tt. At the resonance of (m, 5,) = 
(3,-1/2), 6:^ X n and ^3 = n/2. Then P,(6 = -n/2) shows 
a peak. A similar dip-peak structure of P^ is observed around 
F-wave resonance in Figs. |2] and [3] for our devices. 



P,(0 = -n/2) 



4 sin 62 sin(253 - S2) sin A63 
4 sin^ 62 + sin^ AS^ 



(A24) 



APPENDIX B: CASE OF /tpi^o = 1 



We present calculated results in the case of attractive poten- 
tial and kpRo = 1, where the number of conduction channels 
is unity in the leads. In this case, an interference effect around 
the junction strongly influences the conductance and spin po- 
larization. We discuss the results only for the three-terminal 
device. 

Figures fTOlA) and [TOl'B) show the conductance G± for 
- ±1/2 and spin polarization P,, respectively, when the po- 
tential depth IVol is gradually changed. As seen in Fig.fTOl'A). 
the conductance G± vanishes two times at \V[)\/Ep » 2.7 and 
5.0. The reason why the conductance completely disappears 
cannot be explained by resonant scattering only. This is due 
to an interference effect around the junction area. The value 
of IVol for G+ - and that for G_ = are different from each 
other, in the presence of the SO interaction. As a result, P, = 1 
at G- = and P; = -1 at G+ = 0, as seen in Fig.fJOfB). 

Figure [TOf C) is a grayscale plot of the density of states in 
the junction area, D{E), in the plane of IVol and energy E of 
electron. The band edge of the lowest conduction channel in 
the leads is at E/Ep - 0.616. A virtual bound state of P- 
like state is seen at the Fermi level at the first zero point of 
the conductance (jVoj/iiF ~ 2.7). This indicates a resonant 
scattering via the virtual bound state there. However, we do 
not observe any virtual bound state at the second zero point. 

We observe a perfect spin polarization of P, - +1 in Fig. 
[TOl B). However, the absolute value of the conductance is very 
small when - ±1. Hence, it should be difficult to apply 
this situation to a spin filter. 
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